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NAIE> What are Dynamic Frames 
Navigation and Ancillary Information Facility 


¢ Dynamic reference frames (“dynamic frames' 
for short) have time-dependent orientation. 


— CK- and PCK-based frames are not considered to be 
dynamic frames although they are time-varying. 


— Some examples are given on the next page. 


¢ Dynamic frames are specified using a frames 
kernel (FK). 


¢ The dynamic frames capability is an extension to 
the original SPICE frames subsystem. 


¢ This capability enables SPICE users to conveniently 
use a wide variety of frames that are not "built in" to 
SPICE. 
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Geocentric Solar Ecliptic (GSE) 

Solar Magnetic (SM) 
Spacecraft-centered roll-celestial frame 
Geocentric Solar Magnetospheric (GSM) 
Geomagnetic (MAG) 


— Using constant north centered geomagnetic dipole 
— Using dipole direction defined by time-dependent Euler angles 


Geocentric Solar Equatorial (GSEQ) 

Solar Equatorial frame for any ephemeris object 
Orbital frame for any ephemeris object 

Earth mean equator and equinox of date 

Earth true equator and equinox of date 

Earth mean ecliptic and equinox of date 
Nadir-oriented frame for planetary orbiter 
Radial, tangential, normal (RTN) 


Examples of Dynamic Frames 


Rationale 
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¢ Why should SPICE offer dynamic frames 
implemented using a frames kernel? Instead... 
— auser could build a C-kernel for any frame. 


— SPICE could provide a limited number of "built-in" dynamic 
frames which wouldn't require a frames kernel. 


— users can create their own routines to implement dynamic 
frames. 
¢ Benefits of using a SPICE FK for this purpose 


— Convenience: using a formula rather than a C-kernel avoids C- 
kernel creation, dissemination, storage, and consistency issues. 


— Flexibility: the dynamic frame mechanism enables creation of a 
vast variety of reference frames. 


— Integration: once defined, and once supporting kernels are 
loaded, dynamic frames may be referenced in SPICE API calls. 


— Correct implementation: extensive testing done by NAIF, and 
perhaps further verified by a user, seems less likely to result in an 
error. 
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NAIE> Using a Dynamic Frame 
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¢ Using a dynamic frame in a SPICE-based program is 
straightforward. 
— At program initialization: 
» Load the needed dynamic frame kernel to make the 


frame definition known to SPICE. 


¢ Note: a single dynamic frames kernel likely contains many 
dynamic frame specifications. 


» Load any kernels on which the dynamic frame depends. 


¢ Some dynamic frames require the use of SPK, FK, PCK, 
CK, or SCLK SPICE kernels. 


— Then, in calls to SPICE routines, refer to the dynamic frame 
by name just as you would do with other kinds of frame 
names, such as "J2000." 
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NAIE> Examples of Using Dynamic Frames 
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¢ Find the 6x6 matrix to transform states from the J2000 frame to 
the Geocentric Solar Ecliptic (GSE) frame at the TDB epoch 
given by ET. 


CALL SXFORM( 'J2000', 'GSE', ET1, XFORM ) 


¢ Look up the state of Jupiter relative to the Earth in the GSE 
frame: 


CALL SPKEZR( 'JUPITER', ET1, 'GSE', 
'NONE' , 'EARTH', STATE, LT ) 


¢ You can use dynamic frames any place in SPICE where a frame 
name is used. 


— However, some restrictions apply to use of dynamic frames in 
SPICE kernels (see “Limitations” in the backup slides). 
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Terminology - 1 
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¢ "Frame" is an abbreviation for "reference frame." 


¢ A frame can be defined as a set of three mutually 
orthogonal, unit-length vectors. 


— These vectors are called "basis vectors." The lines containing 
the basis vectors are the "axes" of the frame. 


— The basis vectors indicate the "positive" axis directions; we label 
these vectors +X, +Y, and +Z. The negatives of these vectors are 
labeled -X, -Y, and -Z. 


— We number the axes as follows: 
X=axis1; Y=axis 2; Z=axis 3 
¢ All of the frames we'll deal with are "right-handed": 
this means +Z is the cross product +X x +Y. 


¢ Areference frame's orientation is defined relative to 
another specified frame: the "base frame." 
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Terminology - 2 
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¢ When we Say that a frame is "time-dependent" or "time- 
varying," we mean: 
— The orientation of the frame is time-dependent. 


— Equivalently, the rotation between the frame and its base frame is time- 
dependent. 


¢ By "evaluating" a frame or "evaluating the orientation of a 
frame," we mean computing the rotation between the frame 
and its base frame. 


— Because of their time-dependent nature, an epoch is required in order to 
evaluate a dynamic frame. 


¢ In the SPICE system, frames are considered to have "centers." 


— The center of a frame is always an ephemeris object, something whose 
location can be specified with an SPK file. 


— Frame centers come into play when light time corrections are used: the 
apparent orientation of a time-dependent frame as seen by an observer is 
affected by the one-way light time between the frame's center and the 
observer. 
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NAIE> Terminology - 3 
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When we Say that a vector is "aligned" with another vector, 
we mean that the angular separation between the two vectors 
is zero. 


We use the terms "defining a frame" and "specifying a frame" 
interchangeably. Both refer to creating a frame definition in a 
frames kernel (FK). 


The notation [theta] n 


indicates a frame rotation of theta radians (unless other units 
are specified) about axis n, where n is one of {1, 2, 3}. This 


transformation rotates vectors by —theta radians about axis n. 


Creating a Frames Kernel 
NAIE> Defining Dynamic Frames 
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¢ To define dynamic frames using a frames kernel, a 
fairly detailed understanding of the SPICE dynamic 
frames capability is required. 


¢ A good understanding of the basic SPICE system 
(in particular, the SPK and frames subsystems) is 
also a prerequisite. 
¢ See the Frames Required Reading document for the 
most detailed documentation available. 
¢ The rest of this tutorial is concerned with: 
— explaining the SPICE dynamic frames capability 
— showing how to create dynamic frames kernels 
» We present many dynamic frame definition examples 
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NAIE> Types of Dynamic Frames 
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¢ NAIF envisages providing two types of dynamic 
frames: 
— parameterized 
— scripted 


¢ As of now only the parametrized type is 
implemented within SPICE; the rest of this tutorial 
focuses on this type. 
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¢ Parameterized dynamic frames are defined using formulas. 
— The code implementing the formulas is built into SPICE. 


— The parameters used in the formulas are specified in a frames kernel. 


¢ Parameterized dynamic frames are grouped into frame 
"families". Each family corresponds to a distinct, parameterized 


geometric formula providing a frame definition. The families are: 
— Two-Vector frames 


— Of-date frames 
— Euler frames 
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Defining Parameterized Dynamic Frames 
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¢ Parameterized dynamic frames are defined using 
"keyword=value" assignments in a frames kernel (FK). 


¢ The following items must be specified in the frame definition: 
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Frame name 
Frame ID code 
» The range 1400000-2000000 is reserved for people outside of the NAIF group 
Class (= 5 for dynamic frames) 
Class ID code ( = frame ID code for dynamic frames) 
Frame center (= name or NAIF ID code for central body) 
Frame definition style (= "PARAMETERIZED’') 


Base frame name ( called “relative’’) 


» The frame definition specifies a rotation from the dynamic frame to the base 
(relative) frame. 


Frame family 
Family-specific assignments 


continued on next page 
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continued from previous page 


— Rotation state 


» Possible states are '"ROTATING' and "INERTIAL’. 
¢ A frame is treated as rotating or inertial for the purpose of velocity transformations. 


» The default dynamic frame rotation state is "ROTATING’. 


¢ For rotating two-vector and Euler frames, the rotation state assignment can be 
omitted from the frame definition. 


¢ For "of-date" frames, the frame definition must either specify the rotation state or 
designate the frame as "frozen" at a specified epoch (but not both). 


— Freeze epoch 


» The presence of this optional assignment in a frames kernel indicates that the 
frame orientation, relative to the base frame, is held constant ("frozen") at the 
specified epoch. 


» Most dynamic frames are not frozen. 
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NAIE> What’s Next? 
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¢ In the next approximately 35 pages we'll provide 
implementation details and examples of frames 
belonging to each of the parameterized dynamic 
frame families: 
— Two-Vector frames 
— "Of-Date" frames 
— Euler frames 
— Product frames 
¢ We'll also explain the optional frame attributes 
which can be assigned to frames belonging to the 
above families: 
— “Frozen” 
— “Inertial” 
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¢ Two-vector frames are defined using two time- 
dependent vectors: the "primary" and "secondary" 
vectors. 
— Each vector may be defined by a variety of geometric means: 


» Position vector 

» Target near point vector 
» Velocity vector 

» Constant vector 


¢ The user associates specified positive or negative 
axes of the two-vector frame with the primary and 
secondary vectors. 
— Two-vector frames are always right-handed and have orthogonal 
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axes, so two non-parallel vectors and associations of axes with these 
vectors suffice to define the orientation of a frame. 


Two-Vector Frame Concepts - 2 
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¢ Primary Vector 


— A specified positive or negative axis of the two-vector frame is 
aligned with this vector. 


» The frame kernel creator assigns to this vector one of the axis 
designations { +X, -X, +Y, -Y, +Z, -Z }. 
— Two degrees of freedom of the frame orientation are removed by 


association of an axis with the primary vector. The third degree 
of freedom is the frame's rotation about the primary vector. 


— Example: a frame's -X axis is aligned with the primary vector: 


X Y 
Primary Vector __ 


Z 
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Two-Vector Frame Concepts - 3 
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¢ Secondary Vector 


— A specified positive or negative axis of the two-vector frame is aligned 
with the component of the secondary vector orthogonal to the primary 
vector. 


» The frame kernel creator associates with this vector one of the 
axis designations { +X, -X, +Y, -Y, +Z, -Z }, where the axis is 
orthogonal to that associated with the primary vector. 

— Example, continued: the frame's +Y axis is associated with the 
secondary vector. The component of the secondary vector 
orthogonal to the primary vector is aligned with the frame's +Y axis. 
The secondary vector thus lies in the frame's X-Y plane. 


x 4 


Primary Vector 


Z 


Secondary Vector 


4 Component of secondary 
- vector orthogonal to primary 
- vector 


X 
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Two-Vector Frame Concepts - 4 
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¢ Secondary Vector, continued 


— Typically the secondary vector itself is not orthogonal to the 
primary vector. 
— However, the secondary vector must be linearly independent 
of the primary vector. 
» Near-degenerate geometry can lead to extreme loss of precision. 
¢ This problem can be difficult to diagnose. 


» SPICE enforces independence using a default angular separation 
tolerance of 1 milliradian. The angular separation of the primary 
and secondary vectors may not differ from 0 or Pi radians by less 


than this tolerance. 
» A frame kernel creator can specify a different tolerance value. 
The frame kernel assignment for this is: 
FRAME <frame_ID>_ANGLE SEP TOL = <tolerance> 


where the tolerance is given in radians. 


— Designers of two-vector frames should ensure that the primary 
and secondary vectors can't become nearly parallel for any 
realistic evaluation epoch. 
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NAIE> Two-Vector Frame Concepts - 5 
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¢ In the next several pages we'll describe methods 
for specifying the primary and secondary vectors 
used in creating a dynamic frame. 
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NAIE> Two-Vector Frame Concepts - 6 
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¢ Using a Position Vector 


— Defined by the position of one ephemeris object relative 
to another. The frame kernel creator specifies: 


» the target 
» the observer 
» the aberration correction 


¢ The vector may optionally be corrected for light time and 
stellar aberration. 


— The epoch at which the position vector is computed is 
supplied via a call to a SPICE API: 


» aS an input to an SPK routine, e.g. SPKEZR, SPKPOS. 


» aS an input to a frame system routine, e.g. SXFORM, 
PXFORM. 


— The reference frame relative to which the vector is 
expressed is not specified by the frame kernel creator. 


» SPICE automatically selects this frame. 
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NAIE> Two-Vector Frame Concepts - 7 
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¢ Using a Target Near Point Vector 


— Defined as the vector from an observer to the nearest point ona 
specified extended target body to that observer. The frame 
kernel creator specifies: 


» the target 
» the observer 
» the aberration correction 


¢ The vector may optionally be corrected for one-way light time 
and stellar aberration. 


¢ When one-way light time correction is used, both the position 
and orientation of the target body are corrected for light time. 


— The extended target body is modeled as a triaxial ellipsoid. 
» Size and shape data are given by a PCK. 
— The epoch is supplied via a SPICE API call. 


— The reference frame relative to which the vector is expressed is 
not specified by the frame kernel creator. 


» SPICE automatically selects this frame. 
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Two-Vector Frame Concepts - 8 
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¢ Using a Velocity Vector 


— Defined by the velocity of a target ephemeris object relative to 
an observing ephemeris object. The frame kernel creator 
specifies: 

» the target 
» the observer 
» the velocity reference frame 
¢ This frame may be distinct from the base frame. 


¢ Different velocity frame choices can lead to radically different 
two-vector frame definitions. 


» the aberration correction 


¢ The velocity vector may optionally be corrected for one-way 
light time and stellar aberration. 


¢ Use of light time correction also implies evaluation of the 
velocity vector's frame at a light time corrected epoch: the 
epoch Is corrected for light time between the velocity frame's 
center and the observer, if the velocity frame is non-inertial. 


— The epoch is supplied via a SPICE API call. 
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NAIE> Two-Vector Frame Concepts - 9 
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¢ Using a Constant Vector 
— The vector is constant in a frame specified by the kernel 
creator. 
» The constant vector's frame may be time-dependent. 
» This frame may be distinct from the base frame. 
— The vector may be specified in a variety of coordinate 
systems. 
» Cartesian 
» Latitudinal 
» Right ascension/declination (RA/DEC) 
— An observer may optionally be associated with a constant 
vector for the purpose of defining aberration corrections. 


» The orientation of the constant vector's frame may 
optionally be corrected for one-way light time between the 
frame's center and the observer: if the frame is non- 
inertial, it is evaluated at a light time corrected epoch. 


ntin n nex 
Dynamic Frames continued on next page 
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NAIE> Two-Vector Frame Concepts - 10 
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continued from previous page 


» A constant vector may optionally be corrected for 
stellar aberration due to motion of observer relative to 
solar system barycenter. 

¢ Stellar aberration can be specified without light time 
correction; the string indicating stellar aberration 
correction alone is ‘'S’ 


— The epoch is supplied via a SPICE API call, as for position 
vectors. 


» If the constant vector's frame is time-dependent, that 
frame is evaluated at this epoch, optionally adjusted 
for light time. 
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Two-Vector Frame Examples - 1 
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Nadir-Oriented Spacecraft-Centered Frame 


Primary vector: spacecraft nadir direction 
Y = Z x X, completing the vector. Aligned with nadir frame's -Z axis. 


right-handed frame. 
Y 


Secondary vector: spacecraft velocity 

relative to S/C center of motion in the 

X J2000 frame. Associated with nadir 
frame's +X axis. 


Normalized component of secondary 
vector orthogonal to primary vector is 
aligned with the nadir frame's +X axis. 


See the next page for the FK specifications for this frame. 
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Two-Vector Frame Examples - 2 
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Nadir-Oriented Spacecraft-Centered Frame: 


Frame kernel specification. 


The -Z axis points from the spacecraft toward the closest point on Mars. 


The component of the inertially referenced spacecraft velocity 
vector orthogonal to Z is aligned with the +X axis. 


The +Y axis is the cross product of the +Z axis and the +X axis. 


\begindata 


FRAME <frame_ name> 

FRAME <frame ID>_NAME 

FRAME <frame ID>_CLASS 

FRAME <frame ID> CLASS ID 

FRAME <frame ID> CENTER 

FRAME <frame ID> RELATIVE 

FRAME <frame ID> DEF STYLE 
FRAME <frame_ ID> FAMILY 

FRAME <frame ID> PRI_AXIS 

FRAME <frame ID> PRI VECTOR DEF 
FRAME <frame ID> PRI OBSERVER 
FRAME <frame ID> PRI_TARGET 
FRAME <frame ID> PRI_ABCORR 
FRAME <frame ID> SEC AXIS 

FRAME <frame ID> SEC VECTOR_DEF 
FRAME <frame ID>_SEC OBSERVER 
FRAME <frame ID>_SEC_ TARGET 
FRAME <frame ID>_SEC_ABCORR 
FRAME <frame ID>_SEC_FRAME 
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Definitions 


<frame ID> 

<frame_ name> 

5 

<frame_ID> 
<orbiter ID> 
'J2000' 
'PARAMETERIZED' 
'TWO-VECTOR' 

'_7! 

"TARGET NEAR POINT’ 
<orbiter ID/name> 
"MARS ' 

"NONE ' 

x! 

"OBSERVER TARGET VELOCITY' 
'MARS' i? 
<orbiter ID/name> 

"NONE ' 

'J2000' 
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Two-Vector Frame Examples - 3 
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Spacecraft "View Frame" 


Secondary vector: 
spacecraft position relative to center of motion. 
Associated with view frame's +Y axis in frame 


kernel. 


X = Y x Z, completing the 
right-handed frame. 
("Out of plane" direction) 


X 


Primary vector: spacecraft 


Y velocity relative to center of 
N : motion in J2000 frame. Aligned 
ormalized component of secondary with view frame's Z-axis Ih 
vector orthogonal to primary vector frame kernel. ("Down track" 
is aligned with the view frame's +Y direction) ; 


axis. ("In plane” direction) 


See the next page for the FK specifications for this frame. 


Dynamic Frames 


29 


Two-Vector Frame Examples - 4 
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Spacecraft "View Frame": Frame kernel specification. 


The +Z axis is aligned with the J2000-referenced velocity of the 
spacecraft relative to Mars. 


The component of the spacecraft position orthogonal to +Z is aligned 
with the +Y axis. 


The +X axis is the cross product of the +Y axis and the +Z axis. 


Definitions 


\begindata 


FRAME <frame_ name> 

FRAME <frame_ID>_NAME 
FRAME <frame ID>_CLASS 
FRAME <frame ID>_CLASS ID 
FRAME <frame_ ID>_CENTER 


<frame_ ID> 
<frame_ name> 
5 

<frame_ ID> 
<orbiter ID> 


FRAME <frame ID> RELATIVE 'J2000' 

FRAME <frame ID>_DEF_STYLE 'PARAMETERIZED ' 

FRAME <frame ID> FAMILY 'TWO-VECTOR' 

FRAME <frame ID> PRI AXIS = 'Z' 

FRAME <frame ID> PRI VECTOR DEF = 'OBSERVER_TARGET VELOCITY' 
FRAME <frame ID> PRI OBSERVER = 'MARS' 

FRAME . <frame_ ID>_ PRI TARGET = <orbiter ID/name> 

FRAME <frame ID> PRI _ABCORR = 'NONE' 

FRAME <frame ID> PRI FRAME = 'J2000' 

FRAME <frame ID>_SEC_ AXIS = ty! 

FRAME <frame ID>_SEC_VECTOR_DEF = 'OBSERVER_TARGET POSITION' 
FRAME <frame ID>_SEC_OBSERVER = 'MARS' 

FRAME <frame ID> SEC TARGET = <orbiter ID/name> 

FRAME <frame ID>_SEC_ABCORR = 'NONE' — 


Dynamic Frames 30 


Two-Vector Frame Examples - 5 
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Geocentric Solar Ecliptic Frame (GSE) 


Primary vector: position of sun relative to earth. 
Aligned with GSE frame's +X axis. 


Z=X x Y, 
completing the 
right-handed frame 


Y 


Normalized component 


; : of secondary vector 
Secondary vector: velocity of sun orthogonal to primary 


relative to earth in J2000 frame. : ; : 

Associated with GSE frame's +Y axis in ee cies with the 

frame kernel. ames axis. 
See the next page for the FK specifications for this frame. 
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Two-Vector Frame Examples - 6 
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Geocentric Solar Ecliptic (GSE) frame: 
+X is parallel to the geometric earth-sun position vector. 


+Y axis is the normalized component of the geometric earth-sun velocity 
vector orthogonal to the GSE +X axis. 


+Z axis is parallel to the cross product of the GSE +X axis 
and the GSE +Y axis. 


\begindata 
FRAME GSE = <frame_ID> Definition 
FRAME <frame ID> NAME = 'GSE' 
FRAME <frame ID>_CLASS = 5 
FRAME <frame ID>_CLASS ID = <frame_ ID> 
FRAME <frame ID> CENTER = 399 
FRAME <frame ID> RELATIVE = 'J2000' 
FRAME <frame ID> DEF STYLE = 'PARAMETERIZED' 
FRAME <frame ID> FAMILY = 'TWO-VECTOR' 
FRAME <frame ID> PRI AXIS = 'X! 
FRAME <frame ID> PRI VECTOR _DEF = 'OBSERVER_TARGET POSITION’ 
FRAME <frame ID> PRI OBSERVER = 'EARTH' 
FRAME <frame ID> PRI TARGET = 'SUN' 
FRAME <frame ID> PRI ABCORR = 'NONE' 
FRAME <frame ID>_SEC_AXIS = 'y' 
FRAME <frame ID> SEC VECTOR DEF = 'OBSERVER_TARGET VELOCITY ' 
FRAME <frame ID> SEC OBSERVER = 'EARTH' 
FRAME <frame ID> SEC TARGET = 'SUN' 
FRAME <frame ID> SEC ABCORR = 'NONE' 
FRAME <frame ID> SEC FRAME = 'J2000' 
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Two-Vector Frame Examples - 7 
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Geocentric Solar Magnetospheric Frame (GSM) 


Secondary vector: North geomagnetic centered 
dipole in IAU_EARTH frame. Associated with 
GSM frame's +Z axis in frame kernel. 


Normalized 
component of 
secondary vector 
orthogonal to 
primary vector is 
aligned with the 
Primary vector: position of sun relative to earth. GSM frame’s +Z 
Aligned with GSM frame's +X axis. axis 


Y=Z x X, 
completing the 
right-handed frame 


See the next page for the FK specifications for this frame. 
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Two-Vector Frame Examples - 8 
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Geocentric Solar Magnetospheric (GSM) frame: 
+X is parallel to the geometric earth-sun position vector. 


+Z axis is normalized component of north centered geomagnetic dipole 
vector orthogonal to GSM +X axis. 


+Y completes the right-handed frame. 


\begindata 
FRAME GSM = <frame_ ID> ee 
FRAME . <frame_ ID>_NAME = 'GSM' Definition 
FRAME  <frame_ ID> | __ CLASS 5 


FRAME <frame ID> CLASS ID = <frame_ ID> 
FRAME <frame ID> CENTER = 399 

FRAME <frame ID> RELATIVE = 'J2000' 

FRAME <frame ID> DEF STYLE = 'PARAMETERIZED' 
FRAME <frame ID> FAMILY = 'TWO-VECTOR' 
FRAME <frame ID> PRI AXIS = 'X! 

FRAME <frame ID> PRI VECTOR_DEF = 'OBSERVER_TARGET POSITION’ 
FRAME <frame ID> PRI OBSERVER = 'EARTH' 

FRAME <frame ID> PRI TARGET = 'SUN' 

FRAME <frame ID> PRI _ABCORR = 'NONE' 

FRAME <frame ID> SEC AXIS = 'Z!' 

FRAME <frame ID>_SEC VECTOR _DEF = 'CONSTANT' 
FRAME <frame ID> SEC FRAME = 'IAU EARTH’ 
FRAME <frame ID> SEC SPEC = 'LATITUDINAL' 
FRAME <frame ID> SEC UNITS = 'DEGREES' 

FRAME <frame ID> SEC LONGITUDE = 288.43 

FRAME <frame ID> SEC LATITUDE = 79.54 
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Two-Vector Frame Examples - 9 
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Spacecraft-Centered Roll-Celestial Frame 


Secondary vector: Lock star direction in J2000 frame, corrected 
for stellar aberration due to spacecraft motion. Associated with 
Roll-Celestial frame's +X axis in frame kernel. 


* 


Normalized 
component of 
secondary vector 
orthogonal to 
primary vector is 
aligned with the SCR 
frame’s +X axis 


Primary vector: position of earth relative to 
spacecraft. Aligned with Roll-Celestial frame's 


+Z axis. Y=Z x X, 


completing the 
right-handed frame 


See the next page for the FK specifications for this frame. ae 
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Two-Vector Frame Examples - 10 
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Spacecraft-centered roll-celestial frame: 
+Z is parallel to the geometric earth-sun position vector. 


+X axis is normalized component of star direction orthogonal to Z axis. The star 
direction is corrected for stellar aberration due to motion of the spacecraft. 


+Y completes the right-handed frame. Definitions 


\begindata 
FRAME <frame_ name> 
FRAME <frame ID>_NAME 
FRAME <frame ID>_CLASS 
FRAME . <frame_ ID> | CLASS ID 
FRAME  <frame_ ID> | CENTER 


<frame_ ID> 
<frame_name> 

5 

<frame_ ID> 
<spacecraft ID> 


FRAME . <frame_ ID> | RELATIVE = 'J2000' 

FRAME | <frame_ ID>_ DEF STYLE = 'PARAMETERIZED' 
FRAME . <frame_ ID>_ FAMILY = 'TWO-VECTOR' 
FRAME . <frame_ ID>_ PRI_AXIS = 'Z' 


FRAME . <frame_ ID>_ PRI VECTOR | DEF 
FRAME - <frame_ ID>_| PRI OBSERVER 


'OBSERVER_TARGET POSITION' 
<spacecraft_ ID/name> 


FRAME <frame ID> SEC RA 

FRAME <frame ID>_ SEC DEC 
FRAME <frame ID> SEC OBSERVER 
FRAME <frame ID>_SEC_ABCORR 


<star right ascension in degrees> 
<star declination in degrees> 
<spacecraft ID/name> 

S 


FRAME <frame ID> PRI TARGET = 'EARTH' 
FRAME <frame ID> PRI ABCORR = 'NONE' 
FRAME <frame ID> SEC AXIS = 'X! 

FRAME <frame ID> SEC VECTOR DEF = 'CONSTANT' 
FRAME <frame ID> SEC FRAME = 'J2000' 
FRAME <frame ID> SEC SPEC = 'RA/DEC' 
FRAME <frame ID> SEC UNITS = 'DEGREES' 
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NAIE> "Of-Date" Frames - 1 
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¢ Of-date frames are associated with user-specified 
bodies and are based on user-selected dynamical 
models. 
— Implementations of these models are built into SPICE. 
¢ The currently supported "“of-date" frame families 
are 
— Mean Equator and Equinox of Date 
— True Equator and Equinox of Date 
— Mean Ecliptic and Equinox of Date 


¢ Currently the Earth is the only supported body for 
of-date dynamic frames. 
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NAIE> "Of-Date" Frames - 2 
Navigation and Ancillary Information Facility 


¢ The supported types of models are: 
— Precession 
— Nutation 
— Mean obliquity 


¢ The of-date frame implementation is intended to be 
flexible... 
— The set of supported bodies can grow over time. 
— The set of supported models can grow over time. 


» SPICE is not forever locked into using a single hard- 
heel implementation, such as the 1976 IAU precession 
model. 
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NAIE> "Of-Date" Frames - 3 
Navigation and Ancillary Information Facility 


¢ Mean Equator and Equinox of Date Family 
— For all reference frames in this family... 
» The frame's relationship to the J2000 frame is given 
by a precession model. 


» The frame kernel creator selects a precession model 
from those built into the SPICE software. 


¢ Currently supported only for the Earth, and only the 
1976 IAU precession model (Lieske model) 


» The frame kernel creator must either specify the 
frame's rotation state or must designate the frame 


"frozen" at a specified "freeze epoch” (but not both). 
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“Of-Date" Frames - 4 
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Earth mean equator and equinox of date frame: 


+Z axis is perpendicular to the mean equator of date and points north. 


+X axis is parallel to the cross product of the +Z axis and 
the north-pointing vector normal to the mean ecliptic of date. 


+Y axis completes the right-handed frame. 


\begindata 


FRAME <frame_ name> 

FRAME <frame_ID>_NAME 

FRAME <frame ID> CLASS 

FRAME <frame ID> CLASS ID 

FRAME <frame ID>_ CENTER 

FRAME <frame ID> RELATIVE 

FRAME <frame ID> DEF STYLE 
FRAME <frame ID> FAMILY 

FRAME <frame ID>_PREC_MODEL 
FRAME <frame ID> ROTATION STATE 
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Definitions 


<frame_ID> 
<frame_name> 
5 

<frame_ ID> 
399 

'32000' 
'PARAMETERIZED' 
"MEAN EQUATOR_AND_EQUINOX_OF DATE' 
"EARTH _IAU _ 1976! 


'ROTATING' ee 


This is currently the 
only allowed value for 
precession model. 
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“Of-Date" Frames - 5 


Navigation and Ancillary Information Facility 


¢ True Equator and Equinox of Date Family 
— For all reference frames in this family... 


» The frame's relationship to the J2000 frame is given 
by a precession model and a nutation model. 


» The frame kernel creator selects models from those 
built into the SPICE software. 
¢ Currently supported only for the Earth 


¢ Currently only the 1976 IAU precession model (Lieske model) 
is allowed. 


¢ Currently only the 1980 IAU nutation model is allowed. 
» The frame kernel creator must either specify the 


frame's rotation state or must designate the frame 
"frozen" at a specified "freeze epoch” (but not both). 
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"Of-Date" Frames - 6 


Navigation and Ancillary Information Facility 


Earth true equator and equinox of date frame: 
+Z axis is perpendicular to the true equator of date and points north. 


+X axis is parallel to the cross product of the +Z axis and 
the north-pointing vector normal to the mean ecliptic of date. 


+Y axis completes the right-handed frame. 
\begindata Definitions 


FRAME <frame name> 


= ts <frame_ID> 
FRAME <frame ID>_NAME 


<frame_name> 


FRAME <frame ID>_CLASS 5 

FRAME <frame ID>_CLASS ID <frame_ ID> 
FRAME <frame_ID>_CENTER 399 

FRAME . <frame_ ID>_ RELATIVE '32000' 

FRAME . <frame_ ID>_ DEF_ STYLE 'PARAMETERIZED' 


FRAME | <frame_ ID>_ FAMILY 

FRAME . <frame_ ID>_ PREC MODEL 
FRAME  <frame_ ID> 1 NUT | MODEL 
FRAME . <frame_ ID> | ROTATION _ STATE 


'TRUE_EQUATOR_AND_EQUINOX_OF DATE' 
"EARTH IAU _ 1976' 


'EARTH IAU_1980' — 
"ROTATING ' *“—— => These are currently the 


only allowed values for 
precession model and 
nutation model. 
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“Of-Date" Frames - 7 


Navigation and Ancillary Information Facility 


Mean Ecliptic and Equinox of Date Family 
— For all reference frames in this family... 


» The frame's relationship to the J2000 frame is given 
by a precession model and an obliquity model. 


» The frame kernel creator selects models from those 
built into the SPICE software. 
¢ Currently supported only for the Earth. 


¢ Currently only the 1976 IAU precession model (Lieske model) 
is allowed. 


¢ Currently only the 1980 IAU mean obliquity model is allowed. 
» The frame kernel creator must either specify the 
frame's rotation state or must designate the frame 
"frozen" at a specified "freeze epoch." 
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"Of-Date" Frames - 8 


Navigation and Ancillary Information Facility 


Earth mean ecliptic and equinox of date frame: 


+Z axis is perpendicular to mean ecliptic of date and points toward 


ecliptic north. 


+X axis is parallel to the cross product of the north-pointing 
vector normal to mean equator of date and the +Z axis. 


+Y axis completes the right-handed frame. 


\begindata 


FRAME <frame_ name> 

FRAME . <frame_ ID>_NAME 

FRAME . <frame_ ID> | __ CLASS 

FRAME . <frame_ ID> | CLASS ID 
FRAME . <frame_ ID> | CENTER 
FRAME . <frame_ ID>_ RELATIVE 
FRAME . <frame_ ID>_ DEF STYLE 
FRAME <frame ID> FAMILY 
FRAME <frame ID>_PREC_ MODEL 
FRAME <frame ID>_OBLIQ MODEL 
FRAME <frame ID>_ROTATION_STATE 


Dynamic Frames 


Definitions 


<frame_ID> 
<frame_ name> 
5 

<frame_ ID> 
399 

'J2000' 
'PARAMETERIZED' 
"MEAN | ECLIPTIC . AND ) EQUINOX | OF | DATE ' 
'EARTH | IAU_ 1976' 

'EARTH IAU_1980' ee 


"ROTATING ' 


These are currently the 
only allowed values for 
precession model and 
obliquity model. 
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Euler Frames - 1 


Navigation and Ancillary Information Facility 


¢ Euler frames are defined by a time-dependent 
rotation relative to a base frame. 


— The rotation from an Euler frame to its base frame is given by three 
Euler angles. 


— Each angle is given by a separate polynomial. 
» The polynomials may have different degrees. 


» The independent variable is a time offset, in TDB seconds, 
from an epoch specified by the frame kernel creator. 
» The units associated with the angles are specified by the 
frame kernel creator. 
» The sequence of rotation axes is specified by the frame kernel 
creator. 
¢ The central axis must differ from the other two. 
¢ The rotation from the Euler frame to the base frame is 


[angle_1]axis 1 [angle_2]axis 2 [angle_3]axis 3 
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Euler Frames - 2 


Navigation and Ancillary Information Facility 


¢ Five examples 
— Dynamic version of earth magnetospheric frame (MAG) 


» Latitude and longitude of the north centered geomagnetic 
dipole are given by polynomials. 


— Spinning spacecraft frame 
» The base frame could be a: 
¢ Built-in inertial frame 
¢ C-kernel frame 
¢ Roll-celestial frame (using lock star) 
¢ Nadir frame 


— Topocentric frames for tracking stations for which crustal plate 
motion is modeled 


» The frame rotation keeps the frame orientation consistent with 
the changing station location. 


continued on next page 
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Euler Frames - 3 


Navigation and Ancillary Information Facility 


— Mean or true body equator and body equinox of date frame, 
where the body is a planet or satellite other than the earth 


» The base frame is an IAU_ <body> frame. 


» The Euler frame "removes" the body's rotation about the 
spin axis. 


— Variation on supported "of date" frame 


» An existing supported "of date" frame is used as the base 
frame. 


» Perturbations to the "of date" frame are expressed using 
Euler angles. 
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Euler Frames - 4 


Navigation and Ancillary Information Facility 


As an example, we construct an Euler frame called IAU_MARS EULER. Frame 
IAU_MARS EULER is mathematically identical to the PCK frame named IAU_MARS. 
The PCK data defining the underlying IAU_MARS frame are: 


BODY499 POLE RA = ( 317.68143 
BODY499 POLE DEC = ( 52.88650 
BODY499_ PM = ( 176.630 


-0.1061 O. ) 
-0.0609 0. 
350.89198226 0. ) 


Relative to the angles used to define the IAU_ MARS frame, the angles for our 
Euler frame definition are reversed and the signs negated. Angular units are 
degrees. Rate units are degrees/second, unlike the PCK units of degrees/day. 


angle 3 is 90 + RA 
PCK: angle 2 is 90 - Dec 
angle 1 is PM 


\begindata 
FRAME IAU MARS EULER 
FRAME <frame ID>_NAME 
FRAME <frame ID> CLASS 
FRAME <frame ID> CLASS ID 
FRAME <frame ID>_CENTER 
FRAME <frame ID> RELATIVE 
FRAME <frame ID> DEF STYLE 
FRAME <frame ID> FAMILY 
FRAME <frame ID> EPOCH 
FRAME <frame ID>_ AXES 
FRAME <frame ID>_UNITS 
FRAME <frame ID> ANGLE 1 COEFFS 
FRAME <frame ID> ANGLE 2 COEFFS 
FRAME <frame ID> ANGLE 3 COEFFS 
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angle 1 is -90 - RA 
Euler Frame: angle 2 is -90 + Dec 
angle 3 is - PM 


<frame_ ID> 


'TAU MARS EULER' Definition 


Pes 
499 


"J2000' 

'PARAMETERIZED' 

"EULER' 

@2000-JAN-1/12:00:00 

(313 ) 

"DEGREES ' 

( -47.68143 0.33621061170684714E-10 ) 
( -37.1135 -0.19298045478743630E-10 ) 
( -176.630 -0.40612497946759260E-02 ) 
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Product Frames - 1 


Navigation and Ancillary Information Facility 


¢ Product frames are defined by a product of one or 
more frame transformations. 

— Each factor in the product is a transformation from one frame 
known to the SPICE system to another such frame. The factor 
transformations may be time-dependent. 

¢ Product frames are a generalization of TK (class 4, 
constant-offset) frames. 


¢ A product frame implementing the transformation 
from a frame BASE to a frame PRODUCT has the 


form 


PRODUCT TO 1 TO 2 TO N 
T =a * 7 a ee 
BASE FROM 1 FROM 2 FROM _N 
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NAIE> Product Frames - 2 
Navigation and Ancillary Information Facility 


¢ In the implementation of the product above, the 
factor transformations are applied in right-to-left 
order. 


¢ The from" and to" frames of a product frame 
definition may be completely arbitrary. The only 
restriction on these frames is that the transformation 
from each from" frame to its corresponding to" 
frame must be computable by SPICE at the time the 
product frame is used. 


¢ Because product frames are parameterized dynamic 
frames, limits on recursion depth for dynamic frames 
imply that while the factors may be dynamic frames, 
they may not be dynamic frames that require a level 
of recursion in order to evaluate their orientation. 
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Product Frames - 3 


Navigation and Ancillary Information Facility 


¢ A product frame is specified in a frame kernel using 
assignments of the following form. The "from" and 
"to" frames must be specified by name: 
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FRAME <frame_ name> 

FRAME <frame ID>_NAME 
FRAME <frame ID> CLASS 
FRAME <frame ID> CLASS ID 
FRAME <frame ID>_CENTER 
FRAME <frame ID> RELATIVE 


FRAME <frame ID>_FAMILY 


FRAME <frame ID>_FROM_FRAMES 


FRAME <frame ID>_TO FRAMES 


<frame_ID> 
<frame_name> 

5 

<frame_ID> 

<center ID or name> 
'<frame_name>' 


"PRODUCT' 


( '<from_frame_1>' 
'<from_frame_ 2>' 


'<from_frame N>' ) 


( '<to_frame 1>' 
'<to_ frame 2>' 


'<to_ frame N>' ) 
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NAIE> Frozen Dynamic Frames - 1 
Navigation and Ancillary Information Facility 


¢ A frozen dynamic frame is a "Snapshot" of a 
dynamic frame at a specified epoch. 
— The frame is frozen relative to the base frame specified by 
the frame kernel creator in the frame kernel definition. 
— The rotation from the frozen frame to the base frame is 
constant. 


— The rotation is not frozen with respect to inertial frames 
unless the base frame is inertial. 
— A frame is designated frozen by the presence of a "freeze 
epoch" specification in the frame definition, for example: 
FRAME <FRAME ID>_FREEZE EPOCH = @1949-DEC-31/22:09:46.861901 
— The freeze epoch is specified using SPICE text kernel rules. 
» The “@” syntax is used. 
» The time system is assumed to be TDB. 
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Frozen Dynamic Frames - 2 


Navigation and Ancillary Information Facility 


Frozen version of Earth mean equator and equinox of date frame: 
+Z axis is perpendicular to the mean equator of date and points north. 


+X axis is parallel to the cross product of the +Z axis and the 
vector normal to the mean ecliptic of date. 


+Y axis completes the right-handed frame. 
\begindata Definitions 


FRAME <frame name> 


> ts <frame_ID> 
FRAME <frame ID>_NAME 


<frame_name> 


FRAME <frame ID> CLASS 5 

FRAME <frame ID>_CLASS ID <frame_ ID> 
FRAME _ <frame _ ID>_CENTER 399 

FRAME . <frame_ ID>_ RELATIVE 'J3J2000' 

FRAME . <frame_ ID>_ DEF_ STYLE 'PARAMETERIZED' 


FRAME <frame ID> FAMILY 
FRAME <frame ID> PREC MODEL 
FRAME <frame ID> FREEZE EPOCH 


"MEAN EQUATOR _AND_EQUINOX_ OF DATE' 
"EARTH IAU 1976' 
@1949-DEC-31/22:09:46.861901 
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Inertial Dynamic Frames - 1 


Navigation and Ancillary Information Facility 


¢ Inertial dynamic frames are specified by setting 
the rotation state to 'INERTIAL' in the rotation 
state assignment: 
FRAME_<FRAME_ID>_ROTATION_ STATE = 'INERTIAL' 


— The ‘INERTIAL state implies the frame is treated as inertial for 
the purpose of velocity transformations. 


— The state transformation between any inertial frame and 
"inertial dynamic frame" has a zero derivative block: the state 
transformation matrix has the form 


Se gia C0) | R ( t) 


where R(t) is a time-dependent rotation. 


continued on next page 
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NAIE> Inertial Dynamic Frames - 2 
Navigation and Ancillary Information Facility 


— In contrast, for any rotating frame R(t), the state 
transformation between any inertial frame and R(t) has a 
corresponding matrix of the form 


R(t) | 0 


dR(t)/dt| R(t) 


— The inertial rotation state: 


» simplifies velocity transformations — velocities are 
transformed by a rotation. 


» may be useful for maintaining consistency with other 
dynamic frame implementations. 


» only makes sense if the "inertial" dynamic frame 
actually rotates very slowly! 
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Inertial Dynamic Frames - 3 


Navigation and Ancillary Information Facility 


Inertial version of Earth true equator and equinox of date frame: 
+Z axis is perpendicular to the true equator of date and points north. 


+X axis is parallel to the cross product of the +Z axis and the 
vector normal to the mean ecliptic of date. 


+Y axis completes the right-handed frame. 
\begindata Definitions 


FRAME <frame name> 


= ts <frame_ID> 
FRAME <frame ID>_NAME 


<frame_name> 


FRAME <frame ID> CLASS 5 

FRAME <frame ID>_CLASS ID <frame_ ID> 
FRAME _ <frame _ ID>_CENTER 399 

FRAME . <frame_ ID>_ RELATIVE 'J3J2000' 

FRAME . <frame_ ID>_ DEF_ STYLE "'PARAMETERIZED' 


FRAME | <frame_ ID>_ FAMILY 

FRAME . <frame_ ID>_ PREC MODEL 
FRAME  <frame_ ID> 1 NUT | MODEL 
FRAME . <frame_ ID> | ROTATION _ STATE 


'TRUE_EQUATOR_AND_EQUINOX_OF DATE' 
"EARTH IAU __ 1976' 


'EARTH IAU_1980' — 
' INERTIAL ' *“—— => These are currently the 


only allowed values for 
precession model and 
nutation model. 
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Backup 


Navigation and Ancillary Information Facility 


¢ Numerical Issues 
¢ Limitations 


Dynamic Frames 
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Numerical Issues - 1 


Navigation and Ancillary Information Facility 


¢ Two-vector frame derivatives may be inaccurate. Let R(t) 
represent a time-dependent rotation: 
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— If R(t) depends on CK data, dR(t)/dt may be inaccurate because CK rates 


frequently have low accuracy. 


If R(t) depends on velocity vectors, then dR(t)/dt depends on acceleration 
determined via numerical differentiation. Typically such derivatives suffer 
loss of accuracy. 


» However, if velocities are "well-behaved," numerically derived 
acceleration can be quite good. Example: GSE frame. 


If R(t) depends on position vectors, the velocities associated with those 
vectors by the SPK system may not be mathematically consistent with the 
positions. This can happen for SPK types with separate polynomials for 
position and velocity, such as types 3, 8, 9, and 14. 

If R(t) depends on aberration-corrected vectors, the associated velocities 


may be inaccurate due to accuracy limitations of the aberration 
corrections applied to velocities by the SPK system. 
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Numerical Issues - 2 


Navigation and Ancillary Information Facility 


¢ Recommendations 


— Avoid using aberration corrections in two-vector frame definitions 
if accurate velocity transformations are required. 


— Be aware of the accuracy of the data on which two-vector frames 
are based. 
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Limitations - 1 
Navigation and Ancillary Information Facility 


¢ Simulated recursion: 


— ANSI Fortran 77 doesn't support recursion, so the SPICE dynamic 
frame system implements limited, simulated recursion. 

» Two levels of recursion are supported for selected SPK and 
Frames System routines. 

— Users must avoid requesting "deeper" recursion than the SPICE 
dynamic frame system can support. 

» When defining dynamic frames: 

¢ Choose J2000 as the base frame for two-vector frames. 

¢ Except for Euler frames, avoid using dynamic frames as base 
frames. 

¢ Try to avoid choosing a dynamic frame as the frame associated 
with a velocity or constant vector. 

» In SPK, CK, or PCK kernels, don't use two-vector frames as 
the base frame relative to which ephemeris or attitude data 
are specified. 

¢ "Of-date" or Euler frames are OK for this purpose. 
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Limitations - 2 
Navigation and Ancillary Information Facility 
¢ Run-time efficiency: 
— Dynamic frame evaluation typically requires more computation 
than is needed for CK- or PCK-based frames. 


» For example, evaluation of a two-vector frame may involve 
several SPK calls. 


» Euler frames are an exception: these are fairly efficient as 
long as they don't have a base frame that requires a lot of 
computation to evaluate. 

— To minimize the performance penalty: 

» use J2000 as the base frame for two-vector frames. 

» use the simplest frames possible for association with velocity 
or constant vectors in two-vector frame definitions. 


¢ Prefer non-dynamic frames to dynamic frames and inertial 
frames to non-inertial frames where there is a choice. 
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